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We present a detailed theoretical analysis of Bragg spectroscopy from a Bose-Einstein condensate at T — OK. 
We demonstrate that within the linear response regime, both a quantum field theory treatment and a mean- 
field Gross-Pitaevskii treatment lead to the same value for the mean evolution of the quasiparticle operators. 
The observable for Bragg spectroscopy experiments, which is the spectral response function of the momentum 
transferred to the condensate, can therefore be calculated in a meanfield formalism. We analyse the behaviour 
of this observable by carrying out numerical simulations in axially symmetric three-dimensional cases and in 
two dimensions. An approximate analytic expression for the observable is obtained and provides a means for 
identifying the relative importance of three broadening and shift mechanisms (meanfield, Doppler, and finite 
pulse duration) in different regimes. We show that the suppression of scattering at small values of q observed by 
Stamper-Kum et al, [Phys. Rev. Lett. 83, 2876 (1999)] is accounted for by the meanfield treatment, and can be 
interpreted in terms of the interference of the u and v quasiparticle amplitudes. We also show that, contrary to 
the assumptions of previous analyses, there is no regime for trapped condensates for which the spectral response 
function and the dynamic structure factor are equivalent. Our numerical calculations can also be performed out- 
side the linear response regime, and show that at large laser intensities a significant decrease in the shift of the 
spectral response function can occur due to depletion of the initial condensate. 

PACS numbers: 03.75 Fi 



I. INTRODUCTION 



In 1999 Ketterle's group at MIT reported a set of experi- 
ments in which condensate properties were measured using 
the technique of Bragg spectroscopy [[j} g|. In those experi- 
ments a low intensity Bragg pulse was used to excite a small 
amount of condensate into a higher momentum state, and 
the Bragg spectrum of the condensate was found by measur- 
ing the momentum transfer for a range of Bragg frequencies 
(ui) and momenta (Tiq). That work established Bragg spec- 
troscopy as a tool capable of measuring condensate proper- 
ties with spectroscopic precision. The theoretical analysis of 
the measurements however, gives rise to a number of issues. 
Ketterle and his colleagues assumed that the spectra gave a di- 
rect measurement of the dynamic structure factor, which is the 
Fourier transform of the density-density correlation function, 
and is familiar as the observable in neutron scattering experi- 
ments in super-fluid helium|Q, Q ph. They also attributed the 
suppression of imparted momentum they observed at low q 
values to correlated pair excitations, and quantum depletion 
of the condensate, and speculated [^j that an accurate descrip- 
tion would require a more complete quantum treatment. The 
purpose of the current paper is to develop a theory of Bragg 
spectroscopy that is valid in the regime of the experiments, 
and to use it to make quantitative calculations and to analyze 
the phenomena that occur in this regime. We also investi- 
gate the relationship between the observable of Bragg spec- 
troscopy (i.e. the momentum transferred to the condensate) 
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and the dynamic structure factor, and show that for a trapped 
condensate there is no regime in which one can simply be ob- 
tained from the other. 

We begin in section II within the framework of many body 
field theory and calculate, in the Bogoliubov approximation, 
the linear response of the condensate to an applied Bragg 
pulse. We obtain expressions for the temporal evolution of 
the quasiparticle operators, and show that they have a nonzero 
mean value, i.e. that the quasiparticles are generated as coher- 
ent states. We demonstrate that within a well defined regime 
the mean values of the Bogoliubov operators are identical to 
the amplitudes obtained from a linearized mean-field (Gross- 
Pitaevskii) treatment. The meanfield treatment therefore pro- 
vides a valid description of the experiments in the regime of 
small excitation. 

A number of meanfield theoretical treatments of Bragg 
scattering from condensates have been given. Blakie and Bal- 
lagh have presented a quantitative meanfield description [|7]] 
which confirmed the analysis of the Bragg spectroscopy shift 
given by the MIT group, and provided analytic estimates for 
a number of quantities, including the momentum width of the 
scattered condensate. Stringari and colleagues have also used 
a meanfield description to analyze Bragg spectroscopy ^|. 
and in addition have used the approach to devise schemes for 
measuring quasiparticle amplitudes [ [To|] , and for making spa- 
tially separate condensates interfere Jlip - I n tne current paper 
we use the Gross Pitaevskii formulation of Bragg scattering as 
presented by Blakie and Ballagh |Q] to analyze the behavior 
observed in the Bragg spectroscopy experiments. 

The observable in the experiments is the momentum trans- 
ferred to the condensate, and in section |l| we define a normal- 
ized version of the expectation value of this quantity that we 
call the spectral response function i?(q, ui). For a trapped con- 



densate the momentum transfer can arise from two sources, 
the Bragg beams or the trap itself, which complicates the anal- 
ysis. The MIT group recognized this issue, and applied the 
Bragg laser pulses only for a small fraction of the trap period, 
and then released the trap. However the tradeoff involved in 
minimizing momentum transfer from the trap by using a short 
Bragg pulse significantly compromises the energy selectivity 
of the process. We examine the influence this has on the re- 
lationship between i?(q, oS] and the dynamic structure factor 
5(q, lo), and we show that in the presence of a trap, the evalu- 
ation of S(q, oo) requires i?(q, lo) to be known for all possible 
pulse lengths. 

The central quantity of Bragg spectroscopy is thus the spec- 
tral response function i?(q, oo) and we derive an approximate 
analytic expression for this quantity, incorporating the effects 
of both the meanfield interaction and the finite duration of 
the Bragg pulse, in section IV. In section [v] we use i?(q, lo) 
to characterize our numerical investigations of Bragg spec- 
troscopy, and we consider a wide range of three dimensional 
axially symmetric scenarios, for which we simulate the ex- 
periments using the meanfield (Gross Pitaevskii) equation for 
Bragg scattering [J7J. We verify the validity range of our ap- 
proximate form for i?(q, lo) by comparing it to the full numer- 
ical results, and we identify the regimes in which one or other 
of the mechanisms of: the meanfield interaction, the Doppler 
effect, and the finite pulse duration, dominates the formation 
of the Bragg spectrum. We also show that our approximate 
form for i?(q, 10) will allow a more accurate estimation of the 
momentum width of a condensate than obtained by previous 
analyses. 

Our numerical simulations allow us to calculate the effect 
on Bragg spectroscopy of laser intensities sufficiently large 
that linear response theory no longer holds. We investigate 
cases where the scattered fraction of the condensate is of or- 
der 20%, and show that the depletion of the ground state con- 
densate leads to a significant reduction of the frequency shift, 
which has not been accounted for in previous analyses. We 
also consider the spectral response function from a vortex, us- 
ing two dimensional simulations. Finally in section VI we 



investigate the energy response of a condensate subject to a 
Bragg pulse of sufficiently long duration that individual quasi- 
particle excitations can be resolved. 



II. LOW-INTENSITY BRAGG SCATTERING THEORY 



In this section we calculate the response of the condensate 
to a Bragg pulse within the linear regime, using two distinct 
approaches. In the first of these (section IIA) we use the many 
body field theory formalism in the Bogoliubov approximation, 
to calculate the temporal evolution of the quasiparticle opera- 
tors. In the second approach (section IIB) we use a meanfield 
(Gross Pitaevskii) equation and obtain the amplitudes of the 
linearized response. The two approaches are shown to give 
identical mean results in section IIC. 
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FIG. 1: Bragg spectroscopy of a Bose-Einstein condensate (BEC). 
(a) Two laser beams with wavevectors ki and k2 and frequencies toi 
and u>2 respectively, create a moving optical potential with wavevec- 
tor q = ki — k2 and frequency lo = wi — 102 (see [Q]). (b) Temporal 
behaviour of the Bragg pulse assumed in this paper. 



A. Many-body field theoretic approach 

The many-body Hamiltonian for Nq identical bosons in a 
trap and subject to a time-dependent Bragg pulse can be writ- 
ten 



H = H + Hi(t), 



where Hq is the usual trapped boson Hamiltonian 



(1) 



dr #t 



2m 



Vt(t) 



Uo 
2 



dr \r/ T \I/ T \I/\I/. 
(2) 



Vt{y) is the trapping potential, which we choose to be har- 
monic. The Bragg interaction Hi(t) arises from two over- 
lapping plane wave laser fields, which have equal amplitudes 
but frequency and wave vector differences oo and q respec- 
tively (see Fig. |](a)). The laser fields are treated classically, 
and their interaction with the internal transition of the atoms 
is characterized by a Rabi frequency il(t) (for the combined 
fields at the intensity peaks) and a detuning A which is large 
and essentially the same for both laser fields. In this regime 
the internal structure for the atoms can be eliminated (see 
[Q] for details) so that the field operator ^ refers only to the 
ground internal state, and Hj (t) takes the form 



Hi(t) = J dr& [fiv(t) cos(q • r 
where V(t) = Q 2 (t)/2\A\ (see Fig. |l](a)). 



(3) 
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1. Bogoliubov transformation 



2. Bogoliubov Hamiltonian 



For a highly occupied stationary state (T w OK) the field 
operator can be written in the Bogoliubov approximation as a 
sum of meanfield and operator parts = (\&) + <f>). In this 
paper we mostly consider a ground state, but we also consider 
the case of a central vortex. Following standard treatments 
(e.g. [12, n3[]) we employ a Bogoliubov transformation for 



Applying the transformation in Eq. (Q), to the Hamiltonians 
j) and p) we obtain their Bogoliubov form, i.e. 



Hi 



(14) 



the operator part to write 



*(r,t) 



E 



N Mr)e~ vt + e'^M-"*) x (4) 
kitJe-^tf+btWe^vtir)), 



Hjit) » iff (t) 



where the condensate is represented by the first term, and bi 
and b\ are the quasiparticle destruction and creation opera- 
tors respectively (in an interaction picture with respect to H ). 
The state tpo — \ipo \ exp(iSo), is an eigenstate solution, with 
eigenvalue hfi of the time independent Gross-Pitaevskii equa- 
tion 

t-2 



2m 



V T (r) 



and the functions {v,i, Uj} are the orthogonal quasiparticle ba- 
sis states [14|. These basis states are orthogonal to the con- 



densate mode and are related to the usual (non-orthogonal) 
quasiparticle basis states {uj, Vi] (given below) by projection 
into the subspace orthogonal to the condensate, i.e. 



U; 

* 



where 



= N J dr\ip \ 2 hV(t)cos(q-r- ut) (15) 

i 

xhV(t) cos(q • r — ujt)\xjjo\ + h.c , 

where Eq is the energy of the highly occupied state (see Eq. 
(f75|)). The quasiparticle transformation diagonalizes Hq to 
quadratic order, and we note that the orthogonal basis {Hi, Vi} 
is required for this diagonalization to be valid. In evaluating 
V'o + NoUo\ipo\ 2 ipo , (5) Hf, terms involving products of quasiparticle operators have 

been ignored. This amounts to neglecting Bragg induced scat- 
tering between quasiparticle states, which is of order 1/s/Nq 
smaller than the terms linear in b\ (or bi). Those linear terms 
are of primary interest here, as they describe the scattering be- 
tween the condensate and quasiparticle states which occurs as 
a result of the energy and momentum transfer from the optical 
potential. 

The time dependent exponentials, exp(±iw.jt), which mul- 
tiply the quasiparticle operators in Eq. (Q) account for the free 
evolution due to Hq, and so the Heisenberg equation 



ai\ipo\, 
■a*\ipo\, 



(6) 
(7) 



dr \ipo\ui = - dr \ipo\vi, 



(8) 



d 



(see [0). 

The form of Bogoliubov transformation used in Eq. (Q) ex- 
plicitly includes the phase (So) of ipo, which is convenient in 
cases where ipo is not necessarily a ground state. We note that 
a number of different sign conventions appear in the literature, 
and ours differs from that in Ref. [OJ. We discuss the differ- 
ent conventions in the Appendix. The Bogoliubov-de Gennes 
equation for {m, Vi} are 



Cu. l + N U \iPo\ 2 v t 



i^i •> 

huiiVi, 



(9) 
(10) 



where 



L = 



in- [bie - 

becomes 

= KHf{t) 
This is easily solved to give 

k(t) = bi(0) - fc{t), 
where /3 4 (i) is a c-number, 



, (16) 



huiibi, (17) 
(18) 

(19) 



-|_(V + *VS ) 2 + V T (r) -kfi + 2N Q U \^ \ 

zm. 



with the orthogonality conditions 



dr{mu* - ViV*} 



dr{uiVj — ViUj} 



(11) 



(12) 



(13) 



dr 



0i(t) = -WN Q / dt'V(t')e l 
Jo 

x(u* + i£)cos(q- r-ut')\ip \. (20) 

We see that the Bragg excitation causes the quasiparticle op- 
erators to develop nonzero mean values. Note that this is a 
complete solution of the physics in the linearized regime, from 
which any observable quantities can be computed. 
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3. Initial conditions 

Our derivation so far has been based on a T = OK Bogoli- 
ubov treatment. For this case the initial state is the quasi- 
particle vacuum state (|0)). Had we started with this ini- 
tial condition and considered evolution in the Schrodinger 
picture we would have found that the system evolves as 
|0) — > exp(i6)\{/3i}), where | {/?»}) is a coherent state, i.e. 
bi\{/3i}) = 0M0j}) and 6 is some phase factor (see fllBj]), 

Equation ( |19| ) also provides insight into finite temperature 
cases where most of the atoms are in the condensate. To fully 
treat the finite temperature case, it is necessary to generalize 
the Bogoliubov treatment to account for the thermal depletion 
(e.g. see [|l2[]), requiring the functions ipo, u%> v i an d w i to be 
solved for in a self-consistent manner (e.g. see fll6| , |T7[ |l8[|). 
In that case, Eq. ( |l9| ) for the evolution of % will still apply, 
and thus we see that the initial statistics of 6j(0) are preserved 
and the mean value (pi) is shifted by /3,-. 



where a,j is defined in Eq. (||). Since the quasiparticles form 
a complete set, we may use projection to obtain a set of equa- 
tions for quasiparticle amplitudes that are equivalent to Eq. 
(E3), namely 



c s (t) = -iv(ty 



dr 



O i(jtt-S (r)) 



u ,■(•'"' "■ cos(q • r — ujt)^ 



+v *e' l ^ t - s ^ cos(q • r - " 



(26) 



Because the quasiparticles occupations are all small compared 
to the condensate mode, we can simplify Eq. (|2^) by setting 

*(r,f) = V^Vo|^o(r)| exp(iS' (r) - ifit), yielding 



Ci(t) 



/ dt'V(t')e luja ' / dr{u*+v*) 



x cos(q • r — cut )\ipo\- 



(27) 



This will of course only provide a good solution while the 
quasiparticle occupations all remain small. 



B. Gross-Pitaevskii equation approach 

The Gross-Pitaevskii equation for a condensate of No par- 
ticles subject to a Bragg pulse was derived in [M], i.e. 



h 



U \V\ 2 

+HV(t) cos(q • r - wt)*(r, t) 



--v 2 + y T (r) 

2m 



*(r,*) 
(21) 



where ^(r, t) is the condensate meanfield wavefunction for 
the atoms in their internal ground state, and is normalized ac- 
cording to J dr |^| 2 = N . The condensate wavefunction can 
be expanded in terms of a quasiparticle basis in the form 



(22) 



E 



(ci(t)Uii 



•cf(iKe* 



C. Comparison of approaches 

Direct comparison of the Gross-Pitaevskii and quantum 
field theoretic results is complicated by the fact that they have 
been derived for different basis sets. We note that although the 
Gross-Pitaevskii analysis was carried out using the {ui,Vi} 
basis, it is equally straightforward to use the orthogonal ba- 
sis {Hi, Vi}. Morgan et al. [ p"4| | have shown that (in analysis 
of the Gross-Pitaevskii equation) transforming between these 
bases affects only the ground state population and gives rise 
to no difference in the quasiparticle occupations. 

Connections between the field theory and simple Gross- 
Pitaevskii results based on Eq. ( ]2l| ) can only be expected to 
exist at T = OK, where thermal effects can be ignored 1 . In 
this regime we begin by considering the relevant vacuum ex- 
pectation of the quasiparticle operator from Eq. ( p 7 ^ ) 



where a are the time dependent quasiparticle amplitudes. 
This expansion has been made using the non-orthogonal 
quasiparticle basis, and the ground state has been assumed 
static. This latter assumption will only be valid while ex- 
citation induced by the Bragg pulse remains small. The 
wavefunction decomposition (Eq. (|22|)) transforms the Gross- 
Pitaevskii Eq. (El]) into 



(23) 



_ e -i(So(r)-^) w ^ cog( - q . r _ wt ^ Ti t y (24) 

Morgan et al. [ |l4] ] have shown how to use the orthogonality 
relations of the quasiparticles, Eqs. ( |T2| ) and (|l3]), to project 
out the quasiparticle amplitudes from a condensate wavefunc- 
tion, namely 



(k(t)) = -zv^Vo / dt'V(t')e 
Jo 

x cos(q • r — ujt')\ipo\. 



dr(u*+v*) 



(28) 



This expression emphasizes the coherent nature of the quasi- 
particle states, and is the same as the Gross-Pitaevskii result 
ofEq. (27), with the identification a(i) *-* (bjCt)). We note 
that the apparent difference between equations ( [27| ) and (|28|), 
where the former depends on the matrix elements involving 
{ui, Vi} and the latter on matrix elements involving {Hi, Vi}, 
disappears with the observation that Uj + Vi = u, + Vi (see 
Eqs. (Q) and (^|)). Thus we have verified that the a(t) and 
(bi(t)) are identical. 



-2a* 



dr 



e i(nt-S (.T)) u *q, 



1 We use the word simple here to distinguish our T = OK Gross-Pitaevskii 
equation from more elaborate forms of Gross-Pitaevskii theory which have 



(25) been used to investigate finite temperature effects. For example see jlsj ). 
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where (nj(0)) is the initial occupation, e.g. thermal occu- 
pation at T ^ OK , and is included for generality. In Eq. 
(^) we have chosen to use the {ui, vi \ quasiparticle basis for 
ease of comparison with previous results by other authors (e.g. 

IpLEoLEilEl). 



D. Quasiparticle occupation 

For the physics we consider in this paper the mean occu- 
pation of the i-th quasiparticle level (i.e. (b\bi) or |cj(i)| 2 for 
the many -body or Gross-Pi taevskii methods respectively) is of 
interest, as it is used to make a quasi-homogeneous approxi- 
mation to the spectral response function in Eq. (]58[). Using 
Eq. (jjjjj) we calculate 

To progress, it is necessary to specify the temporal behavior 
of the Bragg pulse. For simplicity we take the pulse shape as 
being square, with V — V p for < t < T p (see Fig. |](b)), 
x cos(q • r - wt')lV>o| + (ft»(0)), (29) then Eq. (§8|) becomes (evaluated at t = T p ) 



(n z (t)) = N / dt' e^'Vit') dr(u*+v*) 



(k(T p )) = -iy/NoV^-W* 



sin((wi - u)T p /2) 



- uj) 



dr«+<) e ^ r |^| 



sin((cjj + uj)T p /2) 
(u)i + uj) 



dr(u*+v*)e-^ r \ip \ 



(30) 



In this expression, the second term (with denominator cOj + uf) 
will typically be significantly smaller than the first term (with 
denominator uji — uj) since uj > 0, so to a good approximation 
we can ignore the second term. Similarly, the quasiparticle 
occupation result (G9h under the same approximation is 



(hi(T p )) 



nV 2 T p N 



dr (u* 



x%- W ,T p ) + (n,(0)), 



in which 



F(uj,T) = 



2sin 2 (wT/2) 

7rT LO 2 



(31) 



(32) 



is a familiar term in time dependent perturbation calculations. 
In particular F(u>i — to, T p ) in Eq. (^Ij) is sharply peaked in 
frequency about uj = u)%, encloses unit area and has a half 
width of a u ~ 2n/T p . In the limit T p — > oo (while V p T p 
remains small) this term can be taken as a (5-function express- 
ing precise energy conservation, so that only quasiparticles of 
energy huj will be excited, i.e. 



lim (h^Tp)) 

l v — *oo 



*V p 2 T p N 



dr(u* 



x6(uj~oji) + (wi(0)>. 



(33) 



III. OBSERVABLE OF BRAGG SPECTROSCOPY 

In this section we outline the experimental procedure of 
Bragg spectroscopy on condensates. We begin by discussing 
the measured observable, which we refer to as the spectral re- 
sponse function. In [111 0] this observable was assumed to be 



a measurement of the dynamic structure factor. We briefly 
review the dynamic structure factor, and discuss why it is in- 
appropriate for these experiments. 



A. Spectral response function 

In the MIT experiments iQ ^] a low intensity Bragg grating 
was used to excite the condensate in-situ for less than a quar- 
ter of a trap period. Immediately following this the trap was 
turned off, the system allowed to ballistically expand, and the 
momentum transfer to the system was inferred by imaging the 
expanded spatial distribution. The experimental signal mea- 
sured is (see rill 0, Eh) 



fl(q,w) 



P(T p )-P(0) 
hq 



where 



(34) 



(35) 



and P is the momentum expectation of the system. We shall 
refer to i?(q, uj) as the spectral response function. For the 
case of a Gross-Pitaevskii wavefunction, the spectral response 
function can be written as 



J dr^*(r,T p ){-ihV)^{r,T p ) 
kq 



(36) 



where it is assumed that the initial condensate has zero mo- 
mentum expectation value. The factors 7 and hq appearing in 
Eqs. ( p4| ) and ([36]) effectively scale out effects of the Bragg 
intensity and duration, the magnitude of momentum transfer, 
and condensate occupation, so that i?(q, u>) can be interpreted 
as a rate of excitation per atom (normalized with respect to 
V 2 ) within the condensate. 
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B. Dynamic structure factor 



quantity 



The dynamic structure factor has played an important role 
in the analysis of inelastic neutron scattering in superfiuid 
4 He. It has facilitated the understanding of collective modes, 
and has enabled measurements of the pair distribution func- 
tion and condensate fraction in that system, as discussed ex- 
tensively in |Q . In those experiments a monochromatic neu- 
tron beam of momentum ftko is directed onto a sample of 4 He 
and the intensity of neutrons scattered to momentum Kk is 
measured, van Hove [^3|| showed that the inelastic scatter- 
ing cross section of thermal neutrons, calculated in the first 
Born approximation, can be directly expressed in terms of the 
quantity 

S(q,w) = i^ e -^|(m|p q |n)| 2 ^-i? m + i?„), 

m.n 

(37) 

which is called the dynamic structure factor (see [||]). In 
Eq. (|3~7j), | to) and E m are the eigenstates and energy lev- 
els of the unperturbed system, Z is the partition function, 
Pq = j dr i&i(r) exp(— iq ■ r)^(r) is the density fluctuation 
operator, and we have taken 



2U2 



fiu> 
q 



h 2 kl 



2m 2to 
ko-k'. 



(38) 
(39) 



Our choice of notation for these quantities is to facilitate com- 
parison between the matrix elements which arise in the dy- 
namic structure factor and Bragg cases. We note that in the 
Bragg context q and uj refer to the wave vector and frequency 
respectively of the optical potential, whereas in a dynamic 
structure factor measurement, Tiq and Tiio are the momentum 
and energy respectively, transferred to the system from the 
scattered probe. 

For the case of a trapped gas Bose condensate, low inten- 
sity off resonant inelastic light scattering [ j24| ] provides a close 
analogy to neutron scattering from 4 He. Csordas et al. [ |20[ ] 
have shown that the cross section for such inelastic light scat- 
tering, with energy and momentum transfer to the photon of 
—Tilo and —Tiq respectively, can be expressed in terms of the 



x (((hi) + l)8(w - Ui) + (hi)5(uj + u)i)) . 



(40) 



This expression generalizes the dynamic structure factor to the 
case of light scattering and applies at finite temperatures in the 
regime of linear response, where the Bogoliubov theory for 
quasiparticles is valid. At T = OK, where thermal depletion 
can be ignored, the dynamic structure factor (^0|) takes the 
form 



S (q,w) = E / dv(u*+v:)e^ r \M 

i 

(T = OK) 



S(u - uJi), 



(41) 



for which a number of approximate forms have been calcu- 
lated (e.g. see[j§[^,|2p). 



The dynamic structure factor and Bragg spectroscopy 

Although the dynamic structure factor and the spectral re- 
sponse function are distinctly different quantities, they do re- 
semble each other strongly. In fact, the matrix elements in the 
T = OK dynamic structure factor in ( ftl) ) resemble those in 
the expression ( |3~3"| ) for the Bragg induced quasiparticle popu- 
lation,^) in the long time limit. It is easy to show, beginning 
from that 



-Sb(q,w) 



{7 > lim 



{hi(T p )} 



(T = OK) (42) 



where (n,-) is given by Eq. ([33]) and 7 is defined in Eq. (|3f^ 
Note that because Eq. ( (42| ) is evaluated at T = OK, we have 
taken the initial occupation, (hi(0)) in Eq. (|33|), to be zero. In 
practice a very long pulse (T p S> 1/ u>t) of sufficiently weak 
intensity (V^Tp -C 1) could be used to excite the quasiparti- 
cles in the regime necessary for Eq. ( p2] ) to hold. 



In [fil ^] the spectral response function i?(q, u>) was assumed to represent a measurement of the dynamic structure factor, 
S(q,uJ). The argument given in those papers was based on the assumption that the Bragg pulse would excite quasiparticles of 
definite momentum foq and the momentum transfer is hence proportional to the rate of quasiparticle excitation. This is in fact not 
so, although we can show that, in a certain sense, the dynamic structure factor and the spectral response function do determine 
each other. We use the result of Brunello et al. who have shown that the momentum imparted can be related to the dynamic 
structure factor according to 

d -^M = -muZZ + 2N K q (£) 2 J *J [S( q , u/) - 5(-q, -.')] (43) 



where the Bragg scattering has been taken to be in the z-direction and P z is the z-component of the momentum expectation. The 
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quantity Z = (X} 7 =i z j) i s me expectation value of the z center of mass coordinate and evolves according to 



dZ = P 1 
dt m 



(44) 



Taking the initial position and momentum expectations to be zero, Eqs. @ and @ can be solved using (|4|) to give the spectral 
response function in terms of the dynamic structure factor as 



iZ(q,w) 



7rT„ 



, cos(w 2 T p ) — cos([o; — ui']T p ) 



This formula can be inverted, but the inversion formulae are different depending on whether lo z is zero or not: 

S'(q, w) — S{— q, — u>) = Jim R(q,iv 7 T p ) for iv z = 0, 

for ui z 7^ 0. 



I R(q.,u,T p )TpdT p 



In the linearized approximation we are using, we can use (^) to show that 

So(q,w) = S(q,w)-5(-q,-w), w > 0, 



(45) 

(46) 
(47) 

(48) 



so that the differencing cancels out finite temperature effects. Beyond the Bogoliubov approximation, there will however be 
residual finite-temperature effects. 



Thus we see it is possible to determine the zero temperature 
dynamic structure factor in a good degree of approximation 
from experiments on a trapped condensate, provided measure- 
ments are performed for a sufficient range of pulse times T p . 
It is clear that this could be a difficult experiment to imple- 
ment, since the spectral response function will drop off faster 
than 1 /T p for large T p , so the measured signal could become 
very small. But in analogy with the results of Sect. VI, we 



would expect that significant information could be obtained 
by measuring for pulse times T p up to about 5 trap periods. 

On the other hand there is no direct connection between the 
spectral response function and the dynamic structure factor 
for any single time T p unless ui z — 0, in which case we must 
take Tp — ► oo. In fact, if we note that significant structure 
in 5(q, u>) is expected on the frequency scale of uj z , then the 
formula (^5j) shows that a smearing over this frequency scale 
is assured by the form of the integrand in (^5|), independently 
of the value of T p . 

One therefore must conclude that the spectral response 
function, not the dynamic structure factor, is the appropriate 
method of analysis for Bragg scattering experiments. How- 
ever, it is in principle possible to determine the zero temper- 
ature dynamic structure factor in a certain degree of approxi- 
mation from these experiments by use of the inversion formula 



IV. QUASI-HOMOGENEOUS APPROXIMATION TO 

R(q,cv) 

In this section we develop an approximation for the spectral 
response function valid for time scales shorter than a quarter 
trap period. Because trap effects are negligible on this time 
scale we employ a quasi-homogeneous approach, based on 



homogeneous quasiparticles weighted by the condensate den- 
sity distribution. This approximation plays an important part 
in the analysis of the numerical results we present in section 

E 



A. Homogeneous spectral response 

The results developed in section II for the trapped conden- 
sates can be applied to a homogeneous system of number den- 
sity n by making the replacements 



u k ( n y k - r /Vv, 

v k (n)e- lk - r /Vv, 



Ui{r) 

/AVMr) _> y/N /V = 



where V is the volume, and 



u k (n) 
u fc (n) 



wf(n) 



2 V CJ fe W u * 

J io% + 2nU Uk/h, 
2^ ' 



(49) 
(50) 
(51) 



(52) 
(53) 

(54) 
(55) 



We have explicitly written these as functions of density n for 
later convenience. Using Eqs. (^9|)-(^T]) the expectation of the 
quasiparticle operators ( |30| ) (i.e. at T = OK) resulting from 
Bragg excitation is found to be 



8 



-iV r 



iujT p /2 



sin((w£ + w)T p /2) 



iuTp /2 f sm{{tJ^{n)-uj)T p /2) 
K(n)-«) 



iqr 



-jq r 



(56) 



Since the homogeneous excitations are plane waves, evaluating the spatial integral in Eq. (56) selects out quasiparticles with 
wave vector k = ±q. The occupations (&£. 6k) of these states are 



kV*T p N 



K(n)+ W9 (n)) 2 J P(( W f(n)Tw),T 1 p ), 



(57) 



where F{u>, T) is defined in Eq. (B2h. 



Since a quasiparticle created by S k carries momentum hk, the 
total momentum transferred to the homogeneous condensate 
is 



(p) = ?iq((^q> - (n-q)) 



(58) 



and the spectral response function, as defined in Eq. (|34j), 
becomes 



#(q,w) = 7 



l(P> 



7 ((«q) - («-q)) 



(59) 
(60) 



F((a;f (n) -w),T p ) - (n)+o;),T p ) 



(61) 



where we have used (u q (n) + v q (n)) 2 — oj q /ui q (n) (from 
Eqs. (p2h and (p3|)) to arrive at the last result. 



B. Quasi-homogeneous spectral response function 

The density distribution of a Thomas-Fermi condensate is 
given by 



(62) 



where N{n)dn is the portion of condensate atoms in the den- 
sity range n — > n + (in, and n p = Ti[j,/Uo, is the peak density 
(see [§). 

To approximate the spectral response function for the in- 
homogeneous case we multiply the portion of the condensate 
at density n by the homogeneous spectral response function 
( |oT| ) for a homogeneous condensate (of density n) and inte- 
grate over all densities present, i.e. 



R QH (q,u) = / dnN(n)R(q,uj) 



(63) 



= ( 



lhN Vl 



dn ■ 



V 4 "p / Jo ^q( n ) 

F(u? (n) - U ,T P ) - F(u?(n)+u,T p ) 



We shall refer to Rqh as the finite time quasi-homogeneous 
approximation to the spectral response function, or simply the 
quasi-homogeneous approximation. Ignoring the finite time 
broadening effects and assuming exact energy conservation 
in Rqh with the replacement F(oj, T) — > S(ui), reduces Eq. 
( |64| ) to the simpler line-shape expression 



8w q N U Q 



1- 



2u q N U 



(65) 
, (66) 



where we have adopted the notation, I q , as used in the original 
derivation Equation ( |66[ ) is also known as the local den- 
sity approximation to the dynamic structure factor (see 
and has been used to analyze experimental data in [[j], |2[ . We 
emphasize that with T p < (1/4) Xt rap the (5-function replace- 
ment is unjustifiable, as we verify with our numerical results 
in the next section. 



V. BRAGG SPECTROSCOPY 

Bragg spectroscopy can broadly be defined as selective ex- 
citation of momentum components in a condensate, by Bragg 
light fields. In this section we consider the spectral response 
function as a Bragg spectroscopic measurement, and using 
numerical simulations of the Gross-Pitaevskii equation (EHi 
and the analytic results of the previous section we identify 
the dominant physical mechanisms governing the spectral re- 
sponse behavior. We investigate the spectral response func- 
tion for a vortex and identify parameter regimes in which a 
clear signature of a vortex is apparent. From the full numeri- 
cal simulations we are also able to assess the effect of higher 
laser intensities on the spectroscopic measurements. 



A. Numerical results for -R(q, lo) 
1. Procedure 



The numerical results we present for J?(q, ui), are found 
(64) by evolving an initial stationary condensate state (typically a 
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ground state) in the presence of the Bragg optical potential, 
using Eq. (|2l|). At the conclusion of this pulse, the spectral 
response is evaluated using Eq. ([36]). This differs slightly 
from the typical procedure in the experiments, where the sys- 
tem is allowed to expand before destructive imaging is used to 
measure the condensate momentum. However we have veri- 
fied numerically (in cylindrically symmetric 3D cases) that 
condensate expansion (after the pulse) does not alter the mo- 
mentum expectation value. For each desired value of q and 
lo, we repeat our procedure of evolving ^ according to (|2l|), 
and calculating i?(q, lo) immediately after the optical pulse 
terminates. 

For axially symmetric situations, the simulations are cal- 
culated in three spatial dimensions with q oriented along the 
z-axis. When the initial state is a vortex, the interesting case 
of scattering in a direction orthogonal to the vortex core (lying 
on the z-axis) would break the symmetry requirement, so for 
these cases 2D simulations with q directed along the y-axis 
are used. For convenience we use computational units of dis- 
tance ro = ^h/2mujT', interaction strength wq — TwjtTq\ 
and time to = lo^ 1 ; where lot is the trapping frequency in 
the direction of scattering. 



2. Parameter regimes 
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We use square pulses of intensity and duration such that 
typicall y les s than 1% of the condensate is excited, except in 



section 



i V E where we investigate the nonlinear response of the 
condensate. As long as the amount of excitation is small, we 
verify that the spectral response function i?(q, lo) is indepen- 
dent of V p . However, the shape of i?(q, lo) is dependent on 
the pulse duration (in accordance with the frequency spread 
about lo associated with the time limited pulse) and on the 
magnitude of q. 

The momentum hq defined by 



hq = ^2mn p U , 



(67) 



(i.e. qo = l/£, where £ is the condensate healing length) 
characterizes the division between regimes of phonon and free 
particle-like quasiparticle character. We note that the experi- 
mental results in [lj and |Q] report measurements of i?(q, lo) 
in the free particle and phonon regime, respectively. 



B. Underlying broadening mechanisms 

We present in Fig. || spectral response functions calculated 
using the Gross-Pitaevskii simulations in three spatial dimen- 
sions. The spectral response functions in Figs. ||(a) and (b) are 
for a spherically symmetric ground state with qo = 3.8/ro. 
The q values (of the Bragg fields) were chosen so that the case 
in Fig. ^(a) is in the phonon regime (q = 2/ro), while Fig. 
|](b) is in the free particle limit (q = 16/ro). For comparison, 
we present in Figs. |2|(c) and (d) spectral response functions 
for a different ground state of greater nonlinearity, for which 
<7o = 8.4/ro. Fig. ^(c) is in the phonon regime (q — 3.76/ro), 



FIG. 2: Spectral response function i?(q, lo) of a 3D condensate, (a), 
(b) Spherical condensate, NoUo = W 4 wo, fJ, = 14.2o;t. (c), (d) 
Oblate condensate with trap asymmetry A = ^/8, NoUo = 5.6 x 
10 wo, M = 70.9u>t- Full numerical solution for R shown as solid 
line; the local density approximation (I q ) dash-dot; and the finite 
time local density approximation (Rqh) dashed. Bragg parameters 
are V p = 0.2lot and T v = OA/ lot- 



whereas Fig. ^|(d) is in the free particle regime (q — 31.6/ro). 
In all cases we compare the Gross-Pitaevskii calculation of R 
with the lines-shape Iq and quasi-homogeneous result Rqh 
from Eqs. ( [66) ) and ( |64[ ) respectively. 



1. Mechanisms 

The two ground states used in calculating the results pre- 
sented in Fig. [| are both in the Thomas-Fermi limit (i.e. they 
satisfy the condition /i ^> lot) 2 - We can see from the fig- 
ure that the local density approximation I q Eq. (66) used 



by previous authors does not always give a good description 
of i?(q, lo), whereas our quasi-homogeneous approximation 
Rqh (q, ui) Eq. (Q) is much more accurate. We have in- 
vestigated the accuracy of Rqh over a wide parameter range 
which has allowed us to evaluate the relative importance of the 
three underlying broadening mechanisms which contribute to 



For our numerical simulations we use Gross-Pitaevskii eigenstates calcu- 
lated from Eq. (ph 
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i?(q, ui). The mechanisms and their contributions are as fol- 
lows: 

i) The shift in the excitation spectrum due to the meanfield 
interaction (depending on the local density). The range 
of densities present in the condensate cause a spread 
in this shift. The frequency width associated with this 
spread is proportional to the chemical potential [i, and 
we shall refer to this as the density width - see section 

rv 

ii) The Doppler effect due to the momentum spread of the 
condensate; the Doppler broadened frequency width is 
(Xp ■ q/m, where a p is the condensate momentum width. 

iii) The frequency spread in the Bragg grating due to the 
finite pulse time; the width arising from this effect is 

~ vr/Tp. 

2. Relative importance of mechanisms 

The relative importance of these mechanisms varies accord- 
ing to the parameter regime. In Table | we compare the esti- 
mated values of these widths for the simulations in Fig. ^| 
We see that for the case of Fig. ||(a) both finite time and den- 





density 
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momentum 


Fig. 
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1(d) 


70.6 
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These figures are indicative of typical Bragg spectroscopy re- 
sults for a small to medium size condensate. The state used 
in Figs. ^|(c) and (d) corresponds to about 10 7 Na atoms in 
a 100 Hz trap, and was chosen to match some of the features 
of the experiments reported in [[l], 0|, e.g. the peak density of 
this state is ~ 4 x 10 14 atoms/cirr and the chemical poten- 
tial is ~ 6.7kHz. The momentum values chosen correspond 
to those used to probe the phonon and free particle regimes in 
those papers (Bragg grating formed by 589 nm beams at 14° 
and 180° respectively). Computational constraints mean we 
cannot match the experimental (prolate) trap geometry, but in 
a similar manner to the measurements made in [[TJ we scat- 
ter along a tightly trapped direction, although from a conden- 
sate in an oblate trap of aspect ratio V8- We note that in the 
phonon probing experiments [|[, scattering was performed 
in the weakly trapped direction for imaging convenience. It 
is worth emphasizing that the mechanisms accounted for in 
our approximate response function Rqh depend only on the 
peak density (i.e. h/i/Uo), the magnitude of q and the Bragg 
pulse duration. The reason is that condensates in the Thomas- 
Fermi regime with the same peak density will have identi- 
cal density distributions in either prolate or oblate traps (see 
Eq. (|62|)). Thus the quasi-homogeneous approximation pre- 
dicts the same spectral response function for both. Scattering 
in a tightly trapped direction will enhance momentum effects 
not accounted for in Rqh, since spatially squeezing the con- 
densate causes the corresponding momentum distribution to 
broaden. However, it is apparent from Table |[ that for the 
case of large condensates this momentum effect is relatively 
small, and thus we expect our result in Figs. |](c) an d (d) to 
give a reasonably accurate description of the MIT experiments 



TABLE I: Simple estimate of widths of component mechanisms of 
the spectral response functions in Fig. ^. All quantities expressed in 
units of lot- 

sity effects are important; hence the quasi-homogeneous re- 
sult which includes both of these is in good agreement with 
i?(q, u>), whereas I q (uS) which accounts for density effect 
alone is quite inaccurate. In Fig. ||(b) the value of q is much 
larger, increasing the momentum width to a point where it 
is comparable to the other mechanisms. Since both I q and 
Rqh fail to account for the condensate momentum, neither 
approximation to the spectral response function is in particu- 
larly good agreement with the numerically calculated R. In 
both Figs. [|(c) and (d) the density width is an order of magni- 
tude larger than both the finite pulse and momentum widths, 
and so in this regime the simple line-shape expression I q is 
generally adequate. 

3. Experimental comparison 

Taking typical experimental parameters, the state used in 
Figs. ^|(a) and (b) corresponds to about 2.8 x 10 5 Na atoms 
in a 50 Hz trap, with the Bragg grating formed from 589 nm 
laser beams intersecting at an angle of 5° in (a) and 39° in (b). 



C. Spectral response function of a vortex 

In previous work [ p5[ | we showed how Bragg scattering 
from a vortex can produce an asymmetric spatially selective 
beam of scattered atoms which provides an in-situ signature 
of a vortex. In Fig. |^ we compare the spectral response func- 
tions from two dimensional ground and vortex states, and in 
Table || we summarize the density, finite time and momentum 
effects for the cases in Fig. ||. In the low momentum trans- 
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(b) ground 
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(b) vortex 
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TABLE II: Simple estimate of widths of component mechanisms of 
the spectral response functions in Fig. ^. All quantities expressed in 
units of lot- 

fer case Fig. [|(a), the density width is the most significant 
component of the response function width, although the in- 
creased momentum distribution of the vortex state relative to 
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FIG. 3: Spectral response functions R(q, w) for 2D ground (solid) 
and m — 1 vortex (dots) states, and Rqh(2D) for 2D ground state 
(dashed). Both states have NqUq = 500«Jo, and the spectra are 
calculated with (a) q = 6/ro and (b) q — 14/ro. Bragg parameters 
are V p = 0.2lu t and T p = 0.8/Vr- 



the ground state is reflected in the width of R. We also see that 
Rqh for the 2D ground state is a reasonable approximation to 
the full numerical calculation of R. 

In the large momentum transfer case Fig. |(b), Doppler ef- 
fects (which scale linearly with q) become important and the 
two peaks in the vortex spectral response appear. These indi- 
cate the presence of flow (or momentum components) running 
parallel and anti-parallel to the scattering direction. Zambelli 
et al. [[| have discussed the Doppler effect in detail, and they 
utilized the impulse approximation for the dynamic structure 
factor for the regime where this mechanism dominates. The 
essence of this approximation is to project the trapped conden- 
sate momentum distribution into frequency space, and they 
have applied this to the case of a vortex (in the non-interacting 
long pulse limit). The impulse approximation always predicts 
a double peaked response from a vortex state corresponding 
to the Doppler resonant frequencies for the parallel and anti- 
parallel momentum components. Our numerical results for 
i?(q, lo) show that at low momenta this vortex signature is ob- 
scured by the density and finite time effects. 



D. Shape characteristics of the spectral response function 

The spectral response functions shown in Figs. ||and|^ con- 
tain a large amount of information, however the major prop- 
erties can be well represented by a few numbers that describe 
the overall shape characteristics of these curves. The charac- 
teristics of the response curve that we focus on are the area 
under the curve, the shift of the mean response function and 
the rms frequency width. These can be expressed in terms 
of the first few frequency moments (mj of R at constant q, 
defined as 

m„ = J dujuj n R{ci,uj). (68) 



m sh ift = m 1 /m - u q ; 
and the rms response width is 



m rms = x/m 2 /mo - (mi/m ) 2 



(69) 



(70) 



These characteristics have also been considered in experi- 
mental measurements [[I] where, however, only two val- 
ues of q were used (corresponding to Bragg gratings formed 
by counter-propagating and approximately co-propagating 
beams respectively) and the speed of sound (equivalent to q ) 
was changed by altering the condensate density. 

Numerically we can consider any value q which can be re- 
solved on our computational grids. Here we choose for sim- 
plicity to directly vary q and consider how the spectral re- 
sponse function changes. The initial state for these calcula- 
tions is identical to that used in Figs. |^(a)-(b) (i.e. a spheri- 
cally symmetric ground state with NqUq = 10 4 woX because 
density, temporal and momentum effects are all important for 
the range of q values we consider. This should be contrasted 
with the initial state in Figs, ^(c)-(d), where the density effect 
dominates over the other two effects at both small and large q. 

Our numerical results lead to the following observations: 

i) In Fig. ^(a) the area under the response curve (i.e. 
mo) is considered. The reduction of this from unity is 
characteristic of suppression of scattering for q < qo 
caused by the interference of the u and v amplitudes 
in the expression for the quasiparticle population (e.g. 
see Eq. (57|)). Apart from the discrepancy at very low 
q (which exists in all moments and which we discuss 
more fully below), both I q and Rqh are in qualitatively 
good agreement with R. The finite pulse duration effect 
on mo is to reduce it uniformly from the I q prediction, 
a feature which Rqh represents well. 

ii) The mean shift (^9|) arises because of the Hartree inter- 
action between the particles excited and those remain- 
ing in the condensate, as indicated by the Bogoliubov 
dispersion relation (Q). The approximations Rqh and 
I q are seen to be in good agreement with the full numer- 
ical calculation over most of the range of q, apart from 
q = (see below). At large values of q, the mean shift 
saturates at the value 4/i/7. 

iii) The width prediction of Rqh is in good agreement with 
the full numerical calculation at low q (2 < qro < 6), 
but differs as q increases and the role of the condensate 
momentum distribution increases in importance. I q is 
in poor agreement for all q, because it ignores temporal 
broadening effects. 

iv) In Fig. |] the sharp features in the full numerical calcu- 
lation of the moments of R at extremely low q arise be- 
cause the (ground state) condensate momentum wave- 
function has significant density in the region where the 
majority of the condensate is being scattered. This 
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FIG. 4: Spectral response function properties (see text) for a 3D con- 
densate with NoUo = 10 4 u>o- R(ci,u) (solid), Rqh{<\,u) (dashed) 
and Iq(uj) (dash-dot), (a) the zeroth moment; (b) the mean response 
shift, with the shift 4/z/7 indicated (dotted); (c) the rms width of 
the frequency response. Bragg parameters are V p = 0.2wt and 
T p = OA/lut. 



causes stimulated scattering of the condensate atoms, 
and will significantly affect the spectral response func- 
tion when the amount of condensate initially present in 
the region we are scattering to is similar (or larger) than 
the amount of condensate excited with the Bragg pulse. 
This effect has been ignored in the derivation of Rqh 
by taking to have a zero initial occupation in Eq. 



Evaluating condensate momentum width from the spectral response 
function 

The MIT group have used the rms-width of the spectral re- 
sponse function to determine the Doppler width At^Dop of 
the condensate, and hence the momentum width cr p (= (m/q) 
Awdop)- This allowed them to make the important observa- 
tion that the condensate coherence length is at least as large 
as the condensate size [|l[] (also see [^]). The Doppler width 
was obtained from the spectral response profile by assuming 
that the Doppler effect adds in quadrature to the density and 
temporal effects. The precise details of the procedure used in 
[jl|] are not given (see Fig. 3 of [[![]), however an expression for 



0.6 



0.4 



CL 



0.2 



qOOOOOO 



□□nO aDDDDDDD 



ooo 



10 



14 







FIG. 5: The true momentum width a p of the condensate (solid line) 
in comparison to the estimated value from Eqs. ( |7^ ) (squares) and 
j7l| ) (circles). Other parameters as in Fig. ^. 



the momentum width of the form 



(A/ 9 ) 2 



(AF( W ,r p ))= 



1 



(71) 



is implied in the text. A difficulty that arises in applying Eq. 
( f7l"| ) is that the rms-frequency spread AF(w, T p ) due to the 
finite pulse length is ill defined. Formally F(oj, T p ) has an un- 
bounded rms-width, though the central peak has a half width 
of ir/T p . In order to allow the best possible result from the 
MIT procedure, we have treated AF as an arbitrary parame- 
ter, and fitted the function a„ in Eq.n7lh to the actual conden- 



, p inEq. _ 

sate momentum width of the ground state we use. The result 
for a range of q values is shown plotted as circles in Fig. |[ 
and was obtained with a best fit value AF w 0.86 x ir/T p . 
We can see that Eq. ( f7l"| ) does not give a particularly good 
estimate of the true momentum width, indicating that the as- 
sumption of quadrature contributions of density, momentum 
and finite time effects to the spectral response function is in- 
appropriate. Agreement improves as q increases, as would be 
expected, because the Doppler effect dominates over the other 
mechanisms in that regime. 

Our approximate form Rqh provides a more accurate 
means of extracting the Doppler width. We first assume that 
the finite time and density effects are well accounted for by 
Rqh, which has an an rms-width ARqh- We then assume 
that the overall width is obtained by adding this in quadrature 
to the Doppler width, which gives the following estimate a p 
for the momentum width 



m rms 



(ARqh) 2 



(72) 



in which no fitting parameter is necessary. Results for a p are 
also shown in Fig. || and clearly give a better estimate for the 
true width a v than Eq. ([71]). 
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FIG. 6: Laser intensity effects on the spectral response function, 
(a) Total scattered condensate fraction for q = 6/ro and (lowest to 
highest) V p — 0.51, 0.85, 1.1, 1.3 wt- (b) Shift of peak response 
frequency m s hift as a function of the peak scattered fraction for q = 
3/ro (diamonds), 6/ro (solid) and 9/ro (crosses). Other parameters: 
spherically symmetric condensate, NqUo = W 4 u>o, go = 3.8/ro 
andTp = 0.8/wt. 



E. Scattering beyond the linear regime 

Spectroscopic experiments require that the number of 
atoms scattered be sufficient for the clear identification of mo- 
mentum transfer to the system. In [|lj the light intensity was 
chosen so that the largest amount of condensate scattered (at 
the Bragg resonance) was about 20%. The validity of the lin- 
ear theories (e.g. section II) must be questioned at such large 
fractional transfers, where excitations can no longer be con- 
sidered noninteracting (e.g. see [|l4[|). The numerical GPE 
calculations do however remain valid in this regime, and can 
be used to understand the changes in the spectral response 
function that occur as the scattered fraction increases. 

In Fig. ^| we present results which extend beyond the linear 
regime. In Fig. IJa) we plot the scattered condensate fraction 
[7 _1 _R(q, oj)/Nq] in a sequence of curves for increasing V p . 
We have chosen to use the scattered condensate fraction rather 
than the spectral response function i?(q, uj), since the numer- 
ical value of the scattered fraction indicates whether the mea- 
surement is outside the linear regime. This is shown clearly 
in this sequence of curves, which would simply be scaled by 
the ratios of if they were in the linear regime. Instead, 
we see that the shapes of the curves change as V p increases 
(and the scattered fraction increases) and in particular the peak 



FIG. 7: Energy response function calculated for a two dimensional 
condensate and a Bragg grating with q — 1/ro and duration (a) T v = 
2tt/ujt and (b) T v = IOtt/wt- In both cases the intensity is chosen 
so that 7" 1 w 0.05JV . 



frequency shifts downwards. In Fig. g(b) we consider that 
dependence of the mean response frequency shift (m s hift) on 
the peak scattered fraction of condensate (i.e. the fraction of 
condensate scattered at the Bragg resonance for the curves) 
for three different q values. A number of features are appar- 
ent in those curves. First, the shift decreases below the linear 
prediction as the scattered fraction increases. Second, the de- 
crease is larger for low momentum transfers ( q < qo) than for 
higher momentum transfers. Third, for large enough momen- 
tum transfer, the shift at a given condensate fraction saturates 

We note that in the MIT experiments a peak excitation frac- 
tion of ~ 0.2 was used [Jl[]. From Fig. ^| we see that even for 
the highest values of momentum transfer, (i.e. q — 6/r or 
q = 9/?*o), the shift will be reduced by of order 15% relative 
to the prediction of the linear theory. 



VI. LONG TIME EXCITATION: ENERGY RESPONSE 

As discussed in section EjJ momentum is not a conserved 
quantity in a trapped condensate, and is not a convenient ob- 
servable for the long time limit. Here we investigate the en- 
ergy transfer to the condensate by Bragg excitation (see also 
[15]), since in the unperturbed system (i.e. in the absence of 
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the Bragg grating) energy is a conserved quantity 3 . A possible 
method of measuring the energy transferred is the technique 
of calorimetry, which has been applied by the Ketterle group 
in a somewhat different context Pq ] . 

We define an energy response function Z(q, ui) analogous 
to the (momentum) spectral response function by the equation 



Z(q, uj) = 7 [E[*(r, T p )] - E[*(r, t = 0)]] , (73) 



where 



E[V] = J dr 



2m 



V 2 + V T (r) 



(74) 



is the energy functional and 7 is defined in Eq. ([35]). Thus 
Z(q, uj) is the rate of energy transfer to the condensate from a 
Bragg grating of frequency uj and wave vector q. The results 
of section II for the quasiparticle evolution can be readily ap- 
plied to evaluating Eq. (73[), Using the T = OK results of 
either the vacuum expectation of the Bogoliubov Hamiltonian 
in Eq. ( p^j ) with bi given by (|l9|), or the linearized Gross- 
Pitaevskii result ( |27| ) in the energy functional ( f74| ) gives 



(75) 



in summation of Eq. ( [77] ) are sufficiently broad in frequency 
space to allow a large number of quasiparticles to respond. 
In Fig. gb) Z(q, uj) is shown for the case of a Bragg pulse 
applied for five trap periods (T p = IQtt/ujt). Here the de- 
tailed structure of the energy response function is revealed, 
with individual (possibly degenerate) quasiparticle peaks be- 
ing visible. We emphasize that current momentum response 
experiments (i.e. R(q, uj)) are essentially limited to at most a 
quarter trap period of excitation. 

A response measurement such as shown in Fig. 0(b) reveals 
a wealth of information about the nature of the condensate ex- 
citation spectrum. The frequencies (cjj) of the quasiparticles 
can be determined by the location of the resonant peaks of 
the response function, and these frequency peaks become nar- 
rower and better defined as T p increases. We note that for a 
given peak corresponding to a Bragg frequency ui p , the area 
under the peak is equal to the matrix element 



duj Z(q, uj) 



TlLJi 



dr«+<)e^ r |Vo| 



(78) 



where Suj (~ [uj p — ir/T p , uj p + tv/T p ]) is the frequency width 
of the pulse and the summation is taken over quasiparticles in 
the energy range Suj. 



where Eq = E[y/Naipo] is the initial (ground state) energy 
and (hi(i)) is the Bragg induced quasiparticle occupation ( |29| ) 
at time t. Substituting Eq. (|75|) into Eq. (73 ) gives 



Z{q,uj) 



7 hiJl 



(76) 



As noted above, because energy is a constant of motion of 
the trapped condensate, we can take T p ^> 2it/ujt (though we 
require V p T p to remain small compared to unity for our linear 
analysis to remain valid) and from Eq. ( pT] ) (with (hi(0)) = 
0), we have 



Z{q, u)=Y^ huj * J dr ( u *i 



uj,T p ). 
(77) 



Thus in long duration (weak) Bragg excitation, the only con- 
tribution to Z(q, uj) will come from quasiparticle states with 
energy approximately matching huj (see Eq. (|33|)). 

In Fig. we present results for the energy response func- 
tion of a 2D condensate for two different durations of Bragg 
excitation. Fig. 0(a) shows the case of Bragg excitation ap- 
plied for a single trap period (T p = 2it/ujt). In this case the 
frequency spread of the Bragg pulse is sufficiently wide that 
only a single broad peak of energy absorption is discernible 
from the condensate i.e. the functions F(ui — uji, T) appearing 



3 We ignore the effect of collisional losses from the condensate which should 
be small for the time scales we consider here. 



VII. DISCUSSION 

In this paper we have given a detailed theoretical analysis of 
the phenomenon of Bragg spectroscopy from a Bose-Einstein 
condensate. We began by deriving analytic expressions for 
the evolution of the quasiparticle operators, which contain all 
possible information about the system in the linear response 
regime. We then demonstrated that at T = OK, the mean 
values of the quasiparticle operators are identical to the quasi- 
particle amplitudes obtained by solving the linearised Gross- 
Pitaevskii equation. Thus, for the purpose of calculating the 
observable of the Bragg spectroscopy experiments (the trans- 
ferred momentum), a meanfield treatment is equivalent to the 
full quantum treatment. Consequently, we based our detailed 
analysis of Bragg spectroscopy on the meanfield equation for 
Bragg scattering presented in a previous paper. The central 
object for the experiments is the spectral response function 
for the momentum transfer, i?(q, uj). We derived the relation- 
ship between i?(q, uj) and the dynamic structure factor and 
showed that, contrary to the assumptions of previous analyses, 
there is no regime in which the two quantities are equivalent 
for trapped condensates. 

The results of our numerical simulations of Bragg spec- 
troscopy, which were carried out in axially symmetric three- 
dimensional cases, or in two dimensions, were characterised 
by the behaviour of i?(q, uj). These full numerical solutions 
are accurate for all values of q that can be resolved by the com- 
putational grid, however the computationally intensive nature 
of the calculation made quantitative comparison of our theory 
with the MIT experiments difficult. The analytic approxima- 
tion for the spectral response function, Rqh (q, uj), provides a 
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means of extending the regime of comparison to large conden- 
sates, and systems without axial symmetry, and we showed 
that it accurately represents i?(q, u>), except in those regimes 
where the momentum width is dominant (a p q/m exceeds 
7t/T p and //) or where stimulated effects can occur (hq < a p ). 
It also provided a means for identifying the relative impor- 
tance of three broadening and shift mechanisms (meanfield, 
Doppler, and finite pulse duration). We have shown that the 
suppression of scattering at small values of q observed by 
Stamper-Kurn et al. [g] is accounted for by the meanfield 
treatment, and can be interpreted in terms of the interference 
of the u and v quasiparticle amplitudes. 

A remaining point to emphasize is that our numerical cal- 
culations allowed us to investigate the regime of large laser 
intensities where the linear response condition is invalid. We 
found that a significant decrease in the shift of the spectral re- 
sponse function can occur due to depletion of the initial con- 
densate. 
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where the quasiparticles modes obey the equations (at T 
OK with No particles in the state) 



Cu[ ±] ± NoUol^v, 



hujM 



(±) 



(80) 



- -h^vl'l. (81) 

and C is as defined in Eq. (pT|). The basis choice of Eqs . 
© and © leads to the vacuum expectation value Eq. (E8h 
having the form 



(&*(*)) 



WNo / dt'V(t')e 



f\ AuJit' 



(82) 



x cos(q • r — wt)|^o| 



The most common form of the Bogoliubov transformation 
(e.g. see Jl2] , |l4||), uses the following form for the expansion 
of the operator part of the field 



(83) 



where the ± indices indicate the relative choice of sign be- 
tween the u and v terms in Eq. (83). These quasiparticles 
modes obey the equations (at T = OK with No particles in the 
-00 state) 



BOGOLIUBOV CONVENTIONS 

Several different forms of the Bogoliubov transformation 
have been used in the theoretical description of inhomoge- 
neous Bose-Einstein condensates. These transformations typ- 
ically differ in their choice of sign between the u and v am- 
plitudes, and by the explicit inclusion of phase. To assist 
comparison of our results to the work of others, we summa- 
rize four different definitions, and show how the form of the 
quasiparticle pop ulation result (Eq. (|5|)) is altered by each 
choice. In section [I A 1 we used the orthogonal quasiparticle 
basis, since the Bogoliubov diagonalization of the many-body 
Hamiltonian must be made with excitations that are orthogo- 
nal to the ground state. However, the quasiparticle population 
result, Eq. (^9|), is insensitive to this and so for brevity we will 
use the non-orthogonal basis states. 

In what follows we denote our choice of quasiparticle am- 
plitudes by the notation {it^~\ v^}. These are closely re- 



lated to the form used in [27| by Fetter (which we denote 
{u\ ,v\ }), but differ by a minus sign in the relative phase 
between Ui and u,. Both these conventions explicitly include 
the condensate phase, and the operator part of the field opera- 
tor is expanded in the form 



[u ( t\v)b i {t)±v^ ) *{v)b\{t) 



(79) 



Lu^ ±N U \i> \ 2 vl ±) = hUiii^, 



where L is defined as 

h 2 



— V 2 + V T {v) - fyi + 2N U \M 
Zm 



(84) 
(85) 



(86) 



In this case the mean value of h is 



(bi(t)) = -iVNo / dt'V{t') 



x / dr (14^*^0 ± Vi^'"^^) cos(q • r — ujt) 



(87) 



(see Eq. (|28|)). When tpo is a ground state (i.e. has constant 
phase), ipo can be taken as real in Eq. (Q), i.e. 

(k(t)} - -i^N~o [ dt'V(t')e iUit ' 



x J dr (u- ± ' ) * ± v^*) cos(q • r — ujt)ipQ. 



(88) 
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